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We consider the range of validity of a Hilbert- trans form approach 
in which the measured magnitude of the modulation-transfer-function 
of an optical fiber is used to compute the fiber's impulse response. It 
is argued that a key "minimum- phase assumption" can fail to be sat- 
isfied in important cases, and a few closely related experimental and 
analytical results are presented. 

I. INTRODUCTION 

Pulse dispersion in an optical fiber transmission line limits its infor- 
mation-carrying capacity by limiting the temporal spacing of input 
pulses that can be resolved at the output. The impulse response g(t), by 
which we mean the output power of a fiber excited by a unit impulse of 
optical power, provides the necessary information concerning the dis- 
tortion of the pulse by modal and material dispersion. For a strictly 
monochromatic pulse source, only modal dispersion contributes to the 
distortion. However, with regard to the corresponding measurement 
problems, it is difficult to obtain sufficiently short (<0.5 nsec) and 
monochromatic (< 10 A) input pulses to accurately study fibers with very 
low modal dispersion. Another difficulty is the lack of availability of 
suitable sources that are tunable over a wide range of wavelengths, in- 
cluding wavelengths longer than 1 fim which are of interest for practical 
fiber systems. 

One can obtain g(t) from the modulation frequency transfer function 
(MTF) G(o>), which is the envelope response of the fiber to an incoherent 
optical signal sinusoidally modulated in amplitude at angular frequency 
a). Personick 1 has shown that to the extent that certain reasonable ap- 
proximations hold, g(t) and C7(oj) are a Fourier transform pair. In prin- 
ciple, the MTF can be determined experimentally. The method employed 
in Refs. 2 and 3 uses a xenon lamp and monochromator as a tunable 
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source that can be sinusoidally modulated to high frequency (>1 GHz) 
by an electrooptic modulator. It is straightforward to measure the 
magnitude |0(oj)| of the transfer function 

G(«) = \G(a)\e u M (1) 

using available components, including an RF spectrum analyzer. How- 
ever, because the fibers must be long (~1 km) in order to obtain good 
measurement precision, the measurement of 0(o>) appears to be formi- 
dable. As / = o)/2ir varies from zero to 1 GHz, 0(a>) varies nearly linearly 
with / from zero to 10 4 7r radians for a 1 km long fiber. However, the 
contribution of 0(a)) to pulse dispersion is due to a nonlinear deviation 
A0(o>), on the order of 2ir, from the much larger linear phase shift 0o(<o); 
i.e., 

0(a>) = o (a>) + A0(a>) (2) 

o (a>) = wL/u (3) 

where L is the fiber length and v is an effective envelope velocity taken 
to be independent of w. Hence, direct measurement of phase distortion 
in the presence of the large frequency-dependent 0o(w) could be subject 
to large error as v varies with temperature or other environmental factors. 
One is therefore led 2,3 to consider methods of mathematical computation 
of 0(a>) from |G(to)| using Hilbert-transform theory. 

The main purpose of this note is to report on results which indicate 
that unfortunately the Hilbert-transform approach described in Refs. 
2 and 3 is, in general, not a helpful one for the particular problem at issue, 
even though early experimental results suggested otherwise.t 

Roughly speaking, it is known that by using a Hilbert-transform 
relation the phase can be obtained from the magnitude of a transfer 
function provided that the transfer function is "minimum phase." A 
standard condition (which is by no means sufficient) for "minimum- 
phase behavior" is that the Laplace transform of the impulse response 
that corresponds to the transfer function have no zeros in the closed right 
half-plane. This is in accord with the observation that for an ideal 
waveguide with constant positive delay tq and transfer function Go = 
e~ iuT °, the phase cannot be determined from a knowledge of the function 
| Go(u) | alone, but a transfer function that has the same magnitude as 
that of the waveguide is Gi(o>) = 1 for which the phase is zero for all 
w.tt 



+ The statement on page 1518 of Ref. 2 concerning the possible lack of "approximate 
minimum phase behavior" was motivated by the results of the joint work described 
here. 

+ * The mathematical reason that Go is not a "minimum-phase" function (even though 
e~ ZT ° has no zeros for Re(2) > 0) is that ln|e _ZT °l fails to satisfy a sufficiently strong growth 
condition in the half-plane Re(z) > 0. See Sections 5.1 and 5.3 for related material. 
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For our purposes, the difference between Go(u>) and Gi(u>) is unim- 
portant, because we are willing to ignore a constant time delay that can 
easily be estimated. More generally, it is reasonable to ask if e IWT0 G(w) 
is a minimum-phase function, where to denotes the linear part of the 
delay. An early impulse response measurement on an actual fiber sug- 
gested 2,3 that this might indeed be the case. However, the further ana- 
lytical and experimental study reported on here shows that nonmini- 
mum-phase behavior is likely to arise, and can arise, in important actual 
cases. Some additional closely related material is also presented. 

Methods for circumventing the difficulties described in this note are 
under study, and it is expected that they will be described in a later 
paper. 

II. SOME ANALYTICAL PROPERTIES OF THE MULTIMODE TRANSFER 

FUNCTION 

In the general case, the transfer function of a fiber can be written in 
the form 

G(o))- C T " e-™ T da(T) (4) 

in which the integrator a (t) is a real- valued monotonically nondecreasing 
function of tJ and T a and Tb, which are fixed by the refractive indices 
of core and cladding, are the smallest and largest modal delays, respec- 
tively. 1 "'' Often a(r) is normalized so that 

da(r) = 1 



I 



For a fiber that can propagate n discrete modes without mode mixing, 
(4) becomes 

G(u) = t dj<? ■*"> (5) 

in which each dj is a positive constant that represents the initial exci- 
tation of the ;th mode. Typically, n > 100. Most of our discussion is 
concerned with the important particular case in which (5) holds, and, 
in order to avoid a lack of continuity of the presentation, proofs of the 
results discussed are given in a separate section. We assume that the r/ 
are ordered so that r\ < ti < . . . < r n . 

In (5), G(o)) is the generalized Fourier transform of a finite train of 



t Thus, roughly speaking, da(r) in (4) can be replaced with b(r)dT in which the function 
6(t) is nonnegative and may contain impulses corresponding to discrete modes. See Refs. 
1 and 2 for the relevant background material. 

t+ We are of course assuming that material dispersion can be neglected. 
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not-necessarily-equally-spaced impulses. Let H(z) denote the corre- 
sponding Laplace transform. That is, let H(z) denote 

for all complex z. Of course G(u) = H{iu) for all w. 

A standard Hilbert-transform method 4 - 2 - 3 for determining the phase 
d(u) of a transfer function (u>), from the function | Q (o>) | , when it is 
possible to do so, is to use the formula* 

= 2co f»° ln|6(y)| 

r 

which, roughly speaking, amounts to a direct application of the Hilbert 
transform 



im-* r*f&« ( 6 ) 

7T Jo y z — OT 



lm\j(iu)}=— 2 o dy (7) 

■k Jo y z — or 



Re[/(oO] 

in which f(z) is any complex- valued function of 2 that satisfies certain 
conditions (such as those described in Section 4.1tt) which include the 
condition that f(z) is analytic for Re(z) > 0. 

The way in which (7) is used to obtain (6) is of course to let f(z) be a 
suitably defined single-valued variant of ln[#(z)], in which #(z) is the 
Laplace transform of the time function whose Fourier transform is G (co). 
That is what gives rise to the well-known requirement that fl{z) be 
zero-free in the closed right half-plane. 

2. 1 The zeros of H(z), and related material 

With regard to the location of the zeros ofH(z), according to Propo- 
sition 2 of Section IV: H(z) ^ for Re (z) > if 

di > t dj (8) 

and if (8) is not satisfied, then given an arbitrarily small positive number 
e, and any set of n numbers ti < t% < . . . < t n , we have H (z) = for some 
z in the closed right half-plane for a choice of the tj such that | r, — tj \ 
< t for all;'. In particular, and roughly speaking, unless (8) is satisfied, 
given any set of n delays, there is a set of delays arbitrarily close to those 



f The integrals in (6) and (7) are to be interpreted as Cauchy principal values. 

^Section 4.1 contains an outline of a proof that (7) holds under certain specific condi- 
tions. The derivation given, for instance in Ref. 4, lacks rigor in that, for example, the point 
s in Ref. 4 is initially assumed to be a point internal to a certain contour, while subsequently 
an expression based on that assumption is evaluated for s on the contour. The main reason 
for including the basically tutorial material of Section 4.1 is that it is used to prove another 
result described in this note. 
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delays such that H(z) is not "minimum phase." Notice that it is not 
claimed that// (z) has a zero in the closed right half-plane whenever (8) 
is violated.* However, Proposition 2 does imply that whenever (8) is not 
satisfied it is incorrect to assert that H(z) ^ for Re (2) > when the 
tj are known only to within some positive tolerance e, no matter how 
small t is. Therefore, H(z ) is zero-free in the closed right half-plane, and 
that property is structurally stable in the sense indicated, if and only if 
(8) is met. This result suggests that it would not be surprising to en- 
counter "nonminimum-phase" behavior with fibers for which the total 
power in a sufficient number of the modes corresponding to the delays 
T2, T3, . . . , t„ is considerable. 

An idealized example in which a somewhat analogous conclusion is 
reached is as follows. Consider an rc-mode fiber without mode mixing 
for which the modal delays are equally spaced by 8 sec, so that t, = t'q 
+ j8 for each positive integer j. Assume that the fractional power into 
the yth mode is given by dj = ce^J for all y, in which 7 and c are constants 
with c > 0. Then 

e (y-z&)n _ 1 

H(z) = ce~ ZT 'o ; - 

and the condition that H(z) 9^ for Re(z) > will be satisfied if and only 
if 7 < 0. t+ Of course 7 < means that modes with larger delays have 
smaller excitation. A similar conclusion is reached for the continuum 
mode-mixing caset in which the integrator a (t) of (4) has a continuous 
derivative that is proportional to e yT . 

While the discussion in the preceding paragraphs suggests t% that there 
are important cases in which (6) cannot be used, it certainly does not rule 
out the possibility that there is some other method for determining 6{u>) 
from I G(co)| (which, for example, might possibly exploit the fact that 
the dj are positive). 

In this connection, consider (5). In order to avoid the necessity of in- 
troducing a function equal to e ,u,ri G(a>), assume throughout the re- 
mainder of this section that t\ = (which of course simply provides a 
normalization 1 '*). 

Suppose, for example, that n = 2 and that di = (1 — di) with < d 1 
< 1 and dt * 1/2. Then | G(«) | = d x 2 + (1 - d{p + 2d 1 (l - d x ) cos (ut 2 ) 



* In fact, we show that that claim would be false. 

tt Since (8) is violated when 7 is negative and sufficiently close to zero, we see that a given 
specific H(z) can be zero-free in the closed right half-plane when (8) is violated. 
I This example was suggested by H. E. Rowe. 

* l Little information is available concerning how to accurately specify G(w) for an actual 
fiber using purely analytical methods. Very small geometrical perturbations can have 
significant effects on the impulse response of graded-index fibers. 5 

+ l It is easy to see that without some such normalization, it is is not possible to determine 
0(u) from I <?(&>) I . 
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Fig. 1 — Measured (deconvolved) impulse responses for (a) a Ge02-Si02 fiber with a < 
aopt; (b) a B 2 3 -Si0 2 fiber with a ~ a opt . 



which clearly is unchanged if di is replaced with (1 — d\). Hence given 
only | G(oj)| for all oj, and that (5) with t\ = holds, it is not possible to 
find a unique 0(o>).t 

With regard to results in the opposite direction, when (5) holds and 
(8) is satisfied, it is true that the phase function 0(a>) can be obtained from 
| G(oj)| by using (6). This is proved in Section 5.3. 

III. EXPERIMENT 

For each of two fibers A and B the impulse response was measured by 
injecting 0.4 nsec pulses (2a width) from a GaAs laser (X = 0.9 ^m) and 
observing the pulse distortions after propagation through the fibers. 
Fiber A was 1 km long and had a graded index Ge02-Si02 core with a 
« 1.9 and a O pt(0.9 /im) » 2.0. Fiber B was 1 km long and had a graded 
index B203-SiC>2 core with a « a opt (0.9 /xm) « 1.8. The measured impulse 
responses are shown in Figs, la and b. These impulse-response functions 
were Fourier transformed to obtain |G(u>)| for each case. Then each 
| G(a>) | together with its phase calculated from piecewise-linear formulas 4 
based on (6) was used to calculate a corresponding impulse response. The 
plots are shown in Figs. 2a and b for fibers A and B, respectively. It may 
be seen that the calculated and measured responses agree quite well for 
fiber B with a suitable normalization and translation to bring them into 



t Analogous examples can be given which hold for all n > 2. For instance, let H(z, y) 
denote the expression for H(z) for the idealized exponential -excitation case mentioned 
above, with c chosen to depend on y so that H (0, 7) = 1. It is not difficult to verify that 
we have \H(iw, y)\ = \H (iu, —y)\ for all u for each 7. 
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Fig. 2 — Impulse responses calculated from the magnitudes of the fourier transforms 
of the measured impulse responses in Figs. 1(a) and (b). 

coincidence. On the other hand, the measured and calculated responses 
for fiber A do not agree.* 

The fact that the functions of Figs, la and 2a appear to be approximate 
mirror images of one another suggests that all of the dominant zeros 
associated with the transfer function of fiber A might lie in the right 
half-plane. In this connection, we note that a general function H(z) of 
the form defined in Section II can have zeros in both half-planes. For 
example, with H(z, y) as defined in a preceding footnote, the product 
H(z, l)H(z, — 1), which can be written in the same form as H(z), has 
zeros in both half planes. 

IV. CONCLUSIONS 

With regard to the overall problem of determining the impulse re- 
sponse g(t), direct measurement of the phase appears to be difficult and 
the general use of the "minimum phase" assumption to calculate the 
phase does not appear to be justified. 

Methods for circumventing the difficulties described are under study, 
and it is expected that they will be described in a later paper. 

V. APPENDIX: PROOFS 

5. 1 Outline of a derivation of a well-known formula 

Let 2 be a complex variable with real and imaginary parts x and y, 
respectively, and let z* denote the complex conjugate of z. Let / be a 

* An early computational error led to a reversal of the sign of the time scale for the phase 
calculations as reported in Refs. 2 and 3. Thus, it was erroneously reported thatg(0 for 
measured and calculated responses matched for fiber A. 
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RADIUS e- 




Fig. 3 — Contour in the (x, y) plane. 



complex-valued function of z defined throughout an open subset S of 
the (x, y)- plane that contains the half-plane x > such that: f(z*) = 
f(z)*, Re[/(iy)] is an even function of y, and Im[f(iy)] is an odd function 
ofy. 



Proposition 1. Suppose that/ is analytic on S, and that \f{z)/z \ 
\z\ -*■ » in the half-plane x > 0. Then 



Oas 



Im[f(iu>)]= — P C 

IT JO 



Mm dy 

y 2 — oj 2 



for each real oj > 0, in which P denotes a Cauchy principal value. 

Outline of a Proof. Assume that the hypotheses of Proposition 1 are 
satisfied. Let co > be given, and choose e > such that e < co and / is 
analytic for \z — iu \ < 2e. With p any positive number such that p > to 
+ e, let C denote the contour shown in Fig. 3 which consists of a semi- 
circular arc of radius p, two semicircular arcs of radius e, and a portion 
of they -axis. 

By Cauchy's integral theorem, 



/(io>) = — i -dz 

2iri Jc z — iu> 
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Therefore, 

Im\f(i*))=f-. f- . f [ Z) dz (10) 

2m Jc (z — io))(z + iu) 

Consider separately the following contributions to the right side of (10): 
1 1 the integral along the y-axis from z = ip to z = -ip excluding the two 
gaps due to the semicircular t arcs, the sum 7 2 of the integrals along the 
c arcs, and the integral I3 along the remaining arc of radius p. 
We find at once that 



ir \_Jo Jw+tJ y z — a' 



h-- 

TT 

Using the fact that the integral of (z — iu>)~ 1 over the upper t arc is iV, 
it follows that J 2 = V2 lm\f(w)] + 5(f) in which 6(e) — as e -* 0. Also, 
^3 = Pip) in which /3(p) -* as p -» «. 
Since Im[/(io;)] = /1 + 7 2 + ^3, we have 

Im[/(i«)] = 2/i + 25(e) + 2j8(p) 

from which it is clear that the limit 

lim lim 2I\ 

p—°° t—Q 

exists and is equal to Im[/(tco)]. This completes the outline of the 
proof.* 

5.2 Result concerning the zeros ofH(z) 

Let H(z) denote" 

for all complex z, in which n > 2, dj > for all ;', and the t, are real 
numbers such that t\ < t 2 < . . . < r n . 

Let t\,t2- ■ • ,t n denote any set of n real numbers with the property 
that ti < * 2 < • • • < *«• 

Consider the condition 

di > t dj (11) 

; = 2 

Proposition 2. We have JY(z) * for Re(z) > if (11) holds. If (11) is 
not satisfied, then for any positive e there is a choice of the tj such that 
I tj — tj\ < e for all j and H(z) = for some z with Re(z) > 0. 



f For further results concerning Hilbert transforms, see, for example, Ref. 6. 
tf For the reader's convenience the definition of H(z) is repeated here. 
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Proof. If (11) holds and z = x + iy with x > 0, 



\H(z)\ = 



g-ZTl 



|"di + L d/c-'^-n)! I > e-*n /rf x - f; «/A > 



Suppose now that 

di < t dj (12) 

7=2 

and let c be given. Let 2 = x + iy with y = ttA. Choose ri = t\ — e, and for 
eachy =2,3, ... ,n, choose t/ such that e , -o , < T ;- T i) = — 1 an( J | T; - — ^ | 
< e. We have 



#(2) = e-"i |di - £ dye-'W-n) 1 



(13) 



in which y — t\ > for) > 2. Using (12) and (13), it is clear that there 
is an x > such that H (z) = 0, which completes the proof. 



5.3 A corollary of proposition 1 

Concerning the function H defined in Section 5.2 consider the fol- 
lowing hypothesis. 

Hypothesis: r\ = and (11) is satisfied. 
We notice that if the hypothesis holds, then, for each real o>, Re[H(io))] 

> and we have H{iu>) = \H(i<a)\e i *^ iu \ in which </>(ia>) denotes the 
principal value of tan -1 \Im[H(ioi)]/Re[H(iw)]\. 

Proposition 3. If the hypothesis holds, then 

0(lw) = — P I — \ f-dy 

■k Jo y l — or 

for a? > in which P denotes a Cauchy principal value. 

Proof: Assume that the hypothesis is satisifed. We shall show that 
Proposition 1 can be used. 
Let £0 be a negative number such that 

di> L dje- x ° T i 

j=2 

We see that there is a positive constant 5 such that Re[tf(z)] > 5 for Re (2) 

> xo- Let 0(2) denote the principal value of tan -1 {Im[/f (2)] /Re [#(2)]} 
for Re(2) > xo. 

For each complex 2 such that 2^0, let p(z) denote the principal value 
of In 2. Thus, p[H(z)] = ln|tf(z)| + itf>(z) forRe(2) > x , and p[H(z)] is 
analytic in 2 for Re(2) > xo. 
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